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We develop a complete analytical description of the time evolution of squeezed states of a charged 
particle under the Fock-Darwin Hamiltonian and a time-dependent electric field. This result gen- 
eralises a relation obtained by Infeld and Plebahski for states of the one-dimensional harmonic 
oscillator. We relate the evolution of a state-vector subjected to squeezing to that of state which 
is not subjected to squeezing and for which the time-evolution under the simple harmonic oscilla- 
tor dynamics is known (e.g. an eigenstate of the Hamiltonian). A corresponding relation is also 
established for the Wigner functions of the states, in view of their utility in the analysis of cold-ion 
experiments. In an appendix, we compute the response functions of the FD Hamiltonian to an 
external electric field, using the same techniques as in the main text. 
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I. INTRODUCTION 

The study of the time-evolution of a non-relativistic charged particle in homogeneous magnetic and electric fields 
has a long history in physics. In a quantum context, the treatment of the problem goes back to Darwin^, who 
considered the evolution of a Gaussian wave-packet in a magnetic field, and Fock^ who obtained the eigenenergies 
and eigenstates of a charged particle in an isotropic harmonic potential, subjected to a magnetic field normal to the 
plane of motion. 

If one takes a particle of charge — e and mass m, moving in the xy plane in an harmonic potencial of frequency loq 
and subjected to a magnetic field B = Be^, the Hamiltonian describing the system is given, in the symmetric gauge 
where A — ^ B x r, hy 
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(1) 



where the operators x,px,y,Py obey the canonical commutation relations. 

This simple problem has applications in the context of the Quantum Hall effect^, where disorder and the Coulomb 
interaction also play a crucial role. Another field for which the study of this Hamiltonian has proved fruitful is that of 
quantum dots, where the simple Hamiltonian given by ([T]) seems to give a good account of the I —Vg curves obtained 
when a gate voltage is applied to the quantum dot^, with corrections due to the assymetry of the confining potential 
and to the Coulomb interactions also playing a role. For some types of quantum dots, such as InAs/GaAs quantum 
dots, the agreement between the theoretical and experimental results seems to hold for magnetic fields up to 15 T'^. 

The study of the evolution of Gaussian wave packets also goes back to the first days of quantum mechanics. This 
study was first undertaken by Schrodinger'', Kennard- and also by Darwii*^ in the context of the harmonic oscillator, 
of a free particle and of a particle in constant electric and magnetic fields. This problem continues to attract attention 
to the present day in many contexts, see the review by Dodono\«^. 

Schrodinger considered the time-evolution of a minimal uncertainty state, i.e. a coherent state of the harmonic 
oscillator, in the terminology of Glauber^. These states have a wide range of applications in quantum optics (see 
e.gjiS,), where they act as the quasi-classical states of the EM field, and in quantum field theory, where they are 
the basis of the phase-space path integral^^. Such states and their derivatives have become important in quantum 
information processing in recent years, in the context of the manipulation of cold atoms in traps. It has become 
possible to reconstruct the Wigner function of a coherent state of the center-of-mass of an harmonically bound ion^-. 

Kennard has on the other hand considered the evolution of a more general wavepacket of the harmonic oscilla- 
tor, what is now known as a squeezed state. Important early contributions are those of Husimi^ and Infeld and 
Plebanski-'^'*'^^. These later authors introduced the so-called squeezing operator and established a relation between 
the evolution of initial states for which the time-evolution is known ('unsqueezed states') and states which are derived 
from such initial states by the application of the squeezing operator ('squeezed states'). The generalisation of this 
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relation to the FD Hamiltonian in an homogeneous time-dependent electric field is the main result of the present 
paper. Stole r^^i^^ proved that squeezed-coherent states are unitarily equivalent to coherent states, thus being minimal 
uncertainty states, but that their minimal character is not preserved by time-evolution, although the uncertainty 
periodically assumes the minimum value compatible with Hcisenberg's uncertainty relation. He also showed that the 
squeezing operator as currently written in terms of quadratures operators is exactly of the form given by Infeld and 
Plebahski. Squeezed states play a proemincnt role in quantum optics, see againiS. Recently, they have also become 
important in quantum information processing both in a quantum optics context and also through the manipulation 
of cold atom a^^i^^'^^i^-'^'^^i^^'^'^'^^i^^'^'^ . See in particular reference^'* for an up-to-date report of the state of quantum 
information processing using cold atoms and photons. 

Besides the paper of Darwin already referred, the time-evolution of states of a non-relativistic particle in homo- 
geneous electric and magnetic fields was also considered by Malkin, Man'ko, Trifonov and Dodono v^^i'^^'^^i'^^'^'^ , who 
considered the dynamics of a particle in an homogeneous electromagnetic field in terms of coherent states, obtaining 
an explicit representation of the Green's function, and studied the invariants of the system; by Lewis and Riesenfeld^, 
who also considered the invariants of such a system; by Kim and Weiner— , who considered the evolution of Gaussian 
wave-packets in a magnetic field, subjected to an isotropic harmonic potential (i.e. the FD Hamiltonian), but also to 
sadle-point potentials, which are relevant for tunnelling problems. 

The structure of this paper is as follows: in the next section, we review the notion of squeezing operator and 
generalise the relation of Infeld and Plebahski to states evolving under the Fock-Darwin Hamiltonian. In section [1111 
with a view to applications in the manipulation of cold atoms, we use the relation obtained to establish a relation 
between the Wigner function of different states and apply it to the special case of squeezed-coherent states. In 
section ITVl we present our conclusions. In appendix we present a derivation of the finite frequency permitivity and 
conductivity of the FD Hamiltonian that uses the same operator methods that are used in the main text but which 
lies somewhat outside of the scope of the main text. Finally, in appendix |B1 we derive the original Infeld-Plebahski 
relation through elementary means. 



II. EVOLUTION OF SQUEEZED STATES UNDER THE FD HAMILTONIAN 

We consider the time-evolution of the state | V't )j that obeys the time-dependent Schrodinger equation ihdt \ ipt) = 
H{t) I V't); where H{t) — Ho + Hi{t), with Hq being given by equation ([T]) and where the interaction Hamiltonean 
Hi{t) of the charge with the external electric field is given, in the Schrodinger picture, by Hi{t) = e{Ex{t)x + Ey{t)y). 

If one now expands the squares and groups the different terms of ([T]), one obtains 

^0 = ^ {Pi +pI) + \ mul [x^ +f) + '^L,, (2) 

where ujl = ^ is the giration frequency and ujj^ =^ luq + with = xpy — ypx being the angular momentum 
component along the z axis. One should note that one can write Hq = Iiq + ^ L^, where ho is the Hamiltonean of 
the isotropic harmonic oscillator with frequency cjfl, and also that [ho, L^] ~ 0. 

Given the state vector | V't )j one can define the corresponding state vector | 0( ) = e"R \ tp^), in the interaction 
representation, such that the two vectors coincide at t = 0. This state vector evolves acording to the Hamiltonian 
H™*'{t) = eR ■^"^ Hi{t) e~ ft which given that ho and in Ho commute, one can also write as 

= Hi{t) e-^^' . (3) 

If one now applies the time-dependent rotation, encoded by to Hi(t), followed by the dynamics of the isotropic 
harmonic oscillator, encoded in ho, one obtains for 

Hr\t) = e {E'^{t) [x cos{ujut) + sm{ujRt)]+ E'y{t) [y cos{ujRt)+py sm{ujRt)]} , (4) 

where 

E'^it) = Exit) cos t) + Ey{t) sin t) , (5) 

E'yit) = -Exit) sin t) + Eyit) cos t) , (6) 

are the components of the electric field in the rotated frame. 
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The wave equation for | (t>t ) can be formally integrated in terms of time-ordered products of the integral of JJ™* it) , 
i.e. \ (pt) =T exp du H"^* (u)^ I "00 )• Now, since the commutator of the operators Hl"*{u) at different times 

is a c-number, one can write the time ordered operator above as 



T exp 



duW^\u) 



exp 



X exp 



duH\''\u) 
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2fi? 



Jo 



(7) 



where the second term on the rhs is a phase factor. 

Collecting the several terms, one obtains for the evolution of | ■i/'f ] 



A) 



exp 



/ ■ 2 



X exp 



\2hmujR Jq 



du / dv sm[uoR{u - v) ] [E'^{u)E'M + E'yi^)E'y(^)] 



2h 



exp 



^ ft-o t ) exp 



(8) 



One now assumes, following Infeld and Plebahski^^, that the initial state | V'o ) related to a certain initial state 
I "00 ) I for which the time-evolution under the isotropic harmonic oscillator Hamiltonian ho is known, by 



■00 ) = exp ( j{P°x + P^y ~ Xqp,^ - YoPy ) j exp (-^r {p ■ r + r ■ p) 



(9) 



where the first operator is a translation operator in phase space, with Xq^Yq, , Py being arbitrary real constants, 
and the second operator is the squeezing operator, with r being a real constant that indicates the degree of squeezing. 

Substituting equation ([51) in ([5]) one can combine the operators e~T^ -/o '^^^l" (") and e'^(^^^+^w^~'''-°^==~^''^!' \ since 
the commutator of their exponents is a c-number. This operation merely generates a phase factor, coming from the 
commutator. One can then commute the resulting operator to the left-hand side, through e~^~'"' e~^'''°*, using 
the time-evolution of x,px,yTPy under and under the time-dependent rotation around the z axis. Combining the 
resulting phase factors and operators, we obtain 



exp 1 / du[Ex{u)xc{u) + Ey[u)yc[u)\] exp [ ^{pl{t)x + Py{t)y - Xcit)po, - ycit)Py) 



X exp 



exp 



j^ht] exp f ^r(p-r + f -p) ) | Vo ) 



(10) 



where Xc{t)^ yc{t), Pxi^) Pyi^) ^^e the classical solutions of the equations of motion for the Fock-Darwin problem 
with initial positions Xq and Yq and initial momenta and Py. These solutions are given by 



Xc{t) 



XqCOS (^'^t 



Yo sin (^^t 



1 



nujjR 
e 



^."cos(^^^t)-P«sm(^t 



du {Ex{u) [sin(aj_|_(t — u)) + sin(ct;_(t — u))] 



2mujR Jo 

+Ey{u) [cos(w+(t — u)) — cos((jj_(t — fi))]} , 



(11) 



yc{t) = 



Xq sin + Yq cos (^^t 



cos{ujRt) 



1 



e 



p2 sin (iy + po cos (iy t) ] siii{ujRt) 

du {—Ex{u) [cos(ti)+(t — u)) — cos{uj-{t — u))] 



2mujR 

+Ey{u) [sin(a;+(t — u)) + sin(ci;_(t — m))]} 



(12) 
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and 



V 2 



sin (^i) cos(wflt) 



e 
2 



XqCOS \^ — t 



Yq sin 



sm{ujiit) 



du {Ex{u) [cos{uj+{t — u)) + cos(cj_(t — u))] 



'Ey[u) [sin(cj+(i — u)) — sin(cj_(t — u))]} , 



(13) 



P° cos 



cos(w_Rt) 



e 
2 



Xq sin ( — 



Yq cos I — t 



sin(a;_Rt) 



du {Ex{u) [sin(ti;+(i — u)) — sin(Lj_(t — u))] 



+Ey{u) [cos(a;+(t — u)) + cos(w_(i — u))]} 



(14) 



One can read the classical dielectric permitivity of the system from equations (|11|12|) and one obtains the same 
results as in appendix [XI which shows that the quantum and classical results are identical, as one would expect for 
a linear system. The classical velocities (t) , Vy (t) can be easily computed by derivation of (|11I12[) with respect to 

time or, using the relations v^{t) = ^ — ^ yc{t)) and Vy{t) = ^ {Py{i) + ^Xc{t)), from equations (fTT|) to 

(fT4)) . One can then read the classical conductivity of the system from the resulting expression and this result again 
coincides with that of appendix El i-c the classical and quantum results are the same. Finally, one can also show, 
after derivation of the velocity expressions with respect to time, that the solution given by (fTTj) and (fT2|) obeys the 
classical equations of motion in two dimensions r = — — {E{t) + r x B) — Wq r, with B = B e^. 

The left most operator in equation (jlOp is again a translation operator in phase phase. If one considers the wave- 
function in the coordinate representation, ip{x,y,t) — {x,y\ ip^), and one applies this translation operator to {x,y |, 

followed by the rotation operator e ss" ^ one obtains the following result 



■0(a;,2/,i) = exp ( -— ( p^(t)a;c(0 +p^?;c(0 + e / du[E^{u) x^u) + Ey{u) ydu)] 



xe^p('-[pUt)x+pl{t)y]] (7^^'• (r-r,(t)) | e"* ^" * (p'^+^p) | > , 



(15) 

with rc{t) being the classical solutions of the equations of motion (|lll) and and TZt being the rotation matrix in 
two dimensions by an angle of LJi,t/2. 

For the wave-function in the momentum representation ■!/)(pa;,py, t) — {px,Py\ipt)j oii6 obtains, applying the 
translation operator to {px,Py |, followed again by the rotation operator. 



-ip{px,Py,t) 



exp 



2h 



pl{t)xc{t) +plyc{t) - e / du[Ex{u)xc{u) + Ey(u)yc{u)] 



X exp 



[p.x,{t)+pyy,{t)] \ (7^^'• (p-p^(t)) I e-^^o'e* 



: r (p-f+f -p) 



■00) , 



(16) 



with p'^{t) being the classical solutions of the equations of motion and ([H)) . 

One can now apply the relation of Infeld and Plebahskii^, relating the evolution of ( a;, y | e^i ^ (p r+r p) | ^ 

tothatof | e~-s''°* | -00 ) and the evolution of ( p^r , | e"* (p'^'+'-p) | Vo ) to that of ( p^^, , | e^-^^"* | ^o), 
to equations (llSp and ()16p . respectively. This relation is derived by elementary means in appendix [Bl 

One obtains for (fT5l) the result 



■4'{x,y,t) = n^^{t) exp(^-^ {pl{i)xc{t)^plyc{t) + ej du [Ea:{u) Xc{u) + Ey{u) ydu)] 



X exp 

X 



-{pl{t)x+pl{t)y) exp 



Mr(t) 



irriLUR sinh(2r) sin(2a>/jt) 



(x - x,{t)f + (y - y,{t))' 



(17) 
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where iirit) — \J e"^'" 003^(0; ^t) + e^^ sin^(a;fl;t) and Tr — -j^ arctan(e^'' tan(aj/jt)). 
For equation (fTBl) . one obtains 



■ip{Px,Py,t) = At-rW '^^P (^PSWa^c(i) +PaycW - e j du [E^,{u) Xc{u) + Ey{u)yc{u)] 



X exp 1 --^L 



where H-r{t) — \J e^*" cos2(w_Ri) + e^^r sin^(cji^t) and t_j. = ^ arctan(e^^'' tan(ci;flt)). The expressions p7)) and 

(fT5)) constitute the main result of this paper and generalise those of Infeld and Plebahski as presented ini^ to states 
evolving under the Fock-Darwin Hamiltonian subjected to a time-dependent electric field in the plane of the system. 



III. THE WIGNER FUNCTION FOR A SQUEEZED STATE 

One can also establish a relation between the Wigner function of a state subjected to squeezing in the presence of 
electric and magnetic fields and the Wigner function of an 'unsqueezed state' evolving under the isotropic harmonic 
oscillator dynamics. 



The Wigner function is defined ai 



,36 



^(r-,P,i) = -^^^ / d''ve-^'P-^^*{r~v/2,t)i;{r + vl2,t). (19) 

The Wigner function gives, after integration with respect to the momentum p or to the coordinate r, the coordinate 
or momentum distribution function of the system and can be considered a quantum generalisation of the Boltzmann 
distribution function. However, the Wigner function is not a bonna-fide distribution since it can assume negative 
values. 

Substituting the result of equation p?)) in (fTO|) . one obtains, after a substitution of variable on v 



(20) 



where u{r,p,t) = iij.it)[p — p'^ (t)] — mujfi sinh(2r) sin(2cjj^^) J"^}^^ and where DQ(r,p,t) is the Wigner function of 
the 'unsqueezed' initial state ^^{x^y) evolving under the isotropic harmonic oscillator dynamics in the absence of 
magnetic or electric fields. 

If the 'unsqueezed' initial state | V'o ) in equation ^ is the vacuum of the isotropic harmonic oscillator, then the 
state I f/'o ) is a squeezed-coherent state, evolving under the Fock-Darwin Hamiltonian. In the context of cold ion 
experiments, one can produce squeezed states of the ion center of mass, either by queenching of the trap frequency 
or parametric amplification^'^, as well as multichromatic excitation of the ion^^iSii^. 

In this case, DQ[r,p,t) — e " Bmi^^ ^ and we obtain for the Wigner function of the squeezed-coherent 
state the result 

D(r,p,t) = ^ C~ "'""f'"'" i^-^At)f-iM^ (P-P'^(t))" + | shM2r) sin(2a.Ht)(p-p'=(t)).(r-r,(t)) ^ ^^^^ 

which represents an assymetric Gaussian, whose shape is preserved by the dynamics of the system, rotating in phase- 
space, with its center determined by the solutions of the classical equations of motion (|llj|14p . Note that for r = 0, i.e. 
for a coherent state, the Wigner function decouples into a product of functions that depend only on the coordinates 
or the momenta. Also, note that for the result (pij) is strictly positive, but this is not the case in general, e.g. if we 
had taken \ tpo) to be an excited-state of the harmonic oscillator. 

A group of still pictures of the x,px section of the Wigner function (PT|) is shown in the figures, for a system 
subjected to a right-handed circularly polarised wave, incident along the zz axis and aligned with the xx axis at i = 0. 
The intensity of the field is Eq — 100 Vm~^, with frequency equal to = 1.42 x 10® Hz. The harmonic frequency of 
the trap is ujq = 7.04 x 10 Hz and the mass of the particle m = 1.50 x 10^^^ Kg, is that of a ®i?e+ ion, with the 



6 



applied magnetic field being B = 6.60 T, which gives a giration frequency ul = 7.06 x 10^ Hz and usji = 7.88 x 10^ 
Hz. Finally, the squeezing parameter r — 0.35. 

The integration of expression (|2ip with respect to p or r yields the coordinate or momentum distribution. The 
wave-packets in coordinate or position space are centered around the classical solutions of the equations of motion, the 
uncertainty in e.g. x,px being given by (Aa;^)t = 2mujR {cos^ {LOjit)er'^^ +sm^ {ujRt)e^^) , {A.pl)t = {cos'^ {LUjit)e'^''+ 

sin^(a-'flt)e~^''), i.e. the uncertainties oscillate with period 2tt/ujii. Their product is given by 

{^x')t {ApDt - ^ (1 + sinh2(2r) sm^i2u;Rt)) > ^ , (22) 

in agreement with Heisenberg's uncertainty relation. Note that the uncertainty in x and Px oscillate in opposition, 
i.e. one increases while the other is decreasing. Also note that, unlike a coherent state, the squeezed-coherent state 
is not a minimum uncertainty state for these two canonical variables, except when t = ■ 

IV. CONCLUSIONS 

In this paper, we have considered the time-evolution of general squeezed states evolving under the Fock-Darwin 
Hamiltonian in an homogeneous time-dependent electric field. We have generalised a relation of Infeld and Plebahski, 
between the time-evolution of states of the harmonic oscillator subjected to squeezing and states not subjected to 
squeezing and for which the time-evolution is known, to states evolving under the FD Hamiltonian in two dimensions. 
A corresponding relation was also established for the Wigner functions of the states. Finally, in an appendix, we 
computed the response functions of the FD Hamiltonian to an external electric field, using the same techniques as in 
the main text. 

Acknowledgements: We acknowledge helpful discussions with A. Ferreira. The authors were supported by FCT 
under Grant No. PTDC/FIS/64404/2006. 



APPENDIX A: CALCULATION OF THE DIELECTRIC PERMITIVITY AND OPTICAL 

CONDUCTIVITY 

We will first indicate how to diagonalise the Hamiltonian ([T]) in the operator representation. The diagonalisation 
of this Hamiltonian and the determination of its eigenfunctions was first obtained by Fock^. 

Introducing the annihilation and creation operators a^, fly, aj^, aj, through x — {h/2mujR)^^^ (clx + a-l.), px — 

i{hmujR/2)^^'^ {cLx — al.), y = {h/2mLJiif'^'^ {ciy + a)y), Py — i {hmijR /2)^^'^ (ciy — ap, one can write such Hamil- 
tonian in the form 

Hq = hujR{al CLx + al (iy + 1) + Lz , (Al) 

with Lz — ifi{a!y ax ~ a\ ^y) being the only operator that is non-diagonal in the annihilation and creation operators 
in the above expression. Introducing the 'circular polarisation' operators a+,a„,aJ^,aL through dx — (a+ + a_), 
Qy — (a+ — a_), a|. = (a^ +al_), a^ = {a\_ — a^), one has Lz = h{al^ a+ — aL a-) and one can write Hq as 

Hq = /ia;-|_a^ a+ + fia;_aL S- + hujR , (A2) 

where lo+ ~ lor + oj^ = t^i? — The energy levels are now given in terms of the occupation numbers of the 
modes +, — hy E = huj+n+ + huj-n^ + Hur. 

We will now compute the dielectric permitivity and optical conductivity of the system in the quantum regime by 
considering the Hamiltonian H(t) of the system interacting with an homogeneous time-dependent electric field, as 
given in section [Hi Expressing the operators x and y in terms of a+, a_, a^_, aL, we have the following expression for 
the interaction Hamiltonian Hi (t) , 

Hi{t) = e{h/2mujRy/^[£+{t)a++£^{t)al+£-{t)a^+£+{t)al], (A3) 

where £+{t) = -^{Ex{t) + iEy{t)) and £^{t) = -^{Ex{t) - iEy{t)). 

If I V't ) is a solution of the time-dependent Schrodinger equation, one defines, as above, the state vector \ (pt) = 
eJi^ot I ^j)^ in the interaction representation, such that the two vectors coincide at t = 0, when the field is turned 
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on. Given that in the interaction representation the annihilation and creation operators contained in Hi{t) evolve in 
time through multiplication by a phase factor e^*"± * , (t) is given by 

+ £_(i)e-*'^-*a_ +£+(i)e*"-*aL] . (A4) 

One can now write, as in sectionlTTl the formal solution \ (pt) — T exp /g du HI"* {u)j | -00 )• The time-ordered 

operator above is then written using identity the second term in this product being a phase factor that can be 
discarded when computing expectation values. 
We can write for | tpt ) the result 

I ) = e-T^ ^{z+{t)al-zlit)a+) + {z.{t)al-z-'_it)a.) | ^ ^ ^^5) 

with z± (t) being given by 

z+{t) = --=^ f du£^{u)e^-^\ (A6) 

= - Jl— f duE+{u)e^'^-\ (A7) 
^/2mnhJii Jo 

and where we have discarded the phase factor referred above. The operator D{z^,z^) — 
g(z+(*)a+-z;(*)a+)+(z-(*)a_-z* (t)a_) ^]^g displacement operator for the annihilation and creation operators, 
i.e. a± D — d± + z±, i)^ aj. D — dj. + Zj^. Therefore, using the above representation of \ ^pt), one can show that 
the average value of any operator A(a+, a^, a_, aL), in the Schrodinger representation, is given by 

{A)t - (Vt \ A{a+, dl,d^,dl)\ ii^t) 

= (Vo |£'M(a+e-^'^+*,a^e*'^+*,a_e-^'^-*,a!_e*"-*)L>| Vo) 
= (^0 |i[(a+ + z+W)e"^"+*,(a^ + <W)e^"+*, 

(a_ + z_(t))e-'"-*, {dl + z* (<))e*--*] | Vo ) • (A8) 

In particular, if ^ = /?^a+ + + j3*_a^ + /?-aL, i.e. for a linear function of the annihilation or creation operators, 
such as x,y,Px,Py, which is the case that will concerns us below, one has that 

6{A)t = f3lz+{t)e-"^+' + (3+zlit)e"^+' + I3*_z^it)e-''^-' + P^z*_{t)e"^-\ (A9) 

where S{A)t = {A)t — {A)^ is the difference between the average value in presence and absence of the applied electric 
field. 

Using this result, one can easily show that the induced polarisation in the system Px{t) — ~e6{x)t, Py{t) = 
—e5{y)t, is given by Pi{t) = /g duxij{t — u) Ej{u), where the permitivity of the system is given by 

Xxxit) = Xyy{t) = ^;^^^isin{uj+t) + siii{uj-t)) , (AlO) 

Xxyit) ^ -Xyxi-t) ^ ^——{C0s{uj+t) - COs{LJ^t)) . (All) 

The first equality between the permitivities follows from rotational invariance around the z axis, the second from 
linear response theory. 

The induced current jx{t) = -e (5( "4 >*, = -e (5( -u^^ )t, is related to P^it) and Py{t), by jx{t) = P^it), jy{t) = 
Py{t) and given that Xij(O) = 0, one easily obtains for the conductivity, defined by ji{t) = duaij{t — u) Ej{u), the 
result fTjj (t) = Xij {t) ■ Hence, 

CTxxit) = (Jyy{t) ^ (w+ C0s(w+t ) + W_ C0s(w_ t) ) , (A12) 

(Jxy{t) = ayx{-t) = (w+ sin(a;+t) - w_ sin(a;_t)) . (A13) 
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Obviously and as above, one can also compute it from the definition oi Vx = ^ {px — i = m {Py + i since 
all the operators involved are linear in the annihilation and creation operators. Note that, since u>+ + a;_ = 2wjj, 
CTaj^CO) = cFyy{Q) = e^/m, a result which agrees with the f-sum rule for a single quantum particle. It is interesting to 
consider this system in two limits, namely B (simple harmonic oscilator) and ujq (particle in a magnetic 

field). In the first case, uj+ = w_ = ojq and one obtains axxit) = cryy{t) = ^cos(ajoi), '^xy{t) = ^yx[—t) = (this 
result is obvious, given the lack of transverse response if i? = 0). In the second case, tj+ — ujl, — 0^ uj^ — 

One has that axx{t) = o-yy{t) = ^ cos(wLt), axyit) = (Tyx{-t) = sin(a;ii). 

Finally, let us consider the case in which a constant electric field is turned on at t = 0. In that case, one obtains at 
large times t ^ oo that the current ji = dij(s 0+) Ej, where ^ij{s) is the Laplace transform of crij{t). Performing 
the integrals, one obtains fTxxi^) — f^yyi^) — 0- The limit s ^ requires a bit of care in the transverse conductivity 

case, since axy{s) = -ayx{s) ^ -5;^^ [ jl+s^ ~ Z?^f^)- obtain axyiO) = -^yx{0) = if t^o 7^ 0. However, we 
obtain axy{0) = —ayx{0) = —e/B in the wq = case (w_ = 0). This result is physically simple to understand if one 
realises that, when a harmonic force is present, a constant electric field merely displaces the force center, whether a 
constant magnetic field is present or not (a shift in the origin of the coordinates merely contributes a constant term 
to the vector potential, that can be simply gauged away). Therefore, one will not observe a response of the velocity 
to the electric field in that case. However, in the absence of an harmonic force, the electric field 'pulls' on the giration 
radius center as if it were a free particle and one does observe a transverse response. 

One should again note that the quantum and classical results obtained for the susceptibilities computed above and 
those obtained from the classical equations of motion (|imi4p are identical and, moreover, that the response to the 
electric field is purely linear. This result follows from the fact that the classical equations of motion and their quantum 
counterparts, the Ehrenfest equations, are linear and can therefore be solved with respect to the field and the initial 
conditions. In this respect, the system behaviour is trivial. However, the equality of results between the classical 
and quantum cases is limited to operators that are linear combinations of the coordinates and momenta. In the 
case of non-linear operators, one can still use the methods discussed in this Appendix to study their time evolution. 
Furthermore, if one is interested in the evolution of wave- functions, as discussed in the main text, one should keep 
the phase factors that were discarded in the computation of average values. 



APPENDIX B: THE INFELD-PLEBANSKI IDENTITY 

We give here an elementary demonstration of the relation of Infeld and Plebahskiii. If Sr = e2k (p ''+'' p) is the 
squeezing operator introduced above, it is easy to show^^ that SrvS^ = re^, SrpS^ = pe~^, i.e. Sr is a scale 
transformation operator that preserves the volume of phase space. Using these identities, one can show that 

^,,|r) = e~^'/^\re~-), (Bl) 

Sr\p) - e'-''/2|pe'-), (B2) 

where d is the space dimension (two in this case). 

We now wish to consider the wave-function of a squeezed state evolving under the isotropic harmonic oscillator ho, 
i.e. the matrix element {x,y \ Sr\ ipo) (the discussion in momentum space is completely analogous). Inserting 

a complete set of position eigenstates, one can write this quantity as 

r + co /• + CO 

ip{x,y,t)^ / / dxQdyoKr{x,y,xo,yo;t)i:o{xo,yo) , (B3) 



00 J —00 



where Kr{x, y, Xq, yo; t) ^ {x,y \ e n^"* Sr \ xq, yo)- Using the identity (jBl|) . one has that 

Kr(x, y, xq, yo; t) = e"'' K{x, y, xq e^"" , yo e^""; t) , (B4) 
where K{x,y,xo,yo',t) = {x,y \ e~Ti'''°* \ xo,yo) is the isotropic harmonic oscillator propagator. Note, however, that 



B4[) is valid for an arbitrary one-particle Hamiltonian. The harmonic oscillator propagator is given by*i 

^\ f imujB. f ,w 2 , 2 , 2 , 2n 2{xxo + yyo)\\ 
K{x, y, xo,yo;t)= exp — — cot cj^t (x +y + Xo + yo) 7^ • B5) 



Now, introducing the scaled variables x — x/^r{t), y = y/fj,r{t), with /ir(i) = \/ e cos^ (wijt) + e^*" sin^ (w^t) , and 
Tr — arctan(e^'" tan(wflt)), one can show, using equation (|B4p . that 

Krix,y,xo,yo;t)^fi;\t)e^^ ^ ^ +^ )k xo,yo;rr] . (B6) 
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Substituting (jB6p in (|B3p . one obtains 

lPix,y,t)=^,-\t)e-''^^ ^■^^H2r)s-n.i2..t)i.^W)^(^_±_^ , (B7) 

which is the Infeld-Plebahski relation used in the main text. Since the expression for the propagator in momentum 
space is entirely analogous to (|B5[) . the steps are identical to those above, except that r is replaced by —r. One 
obtains 



10 



REFERENCES 



14 



Electronic address: "jaime. santos@fisica.uminho.pt" 

C. G. Darwin, "Free Motion in the Wave Mechanics" Proc. Royal Soc. London 117, 258-293 (1927). 
V. Fock, "Bemerkung zur Quantelung des harmonischen Oszillators im Magnetfeld" Z. Phys. 47, 446-448 (1928). 
M. Janssen, O. Viehweger, U. Fastenrath, and J. Hajdu, Introduction to the theory of the integer quantum Hall effect. (VCH, 
1994). 

L. P. Kouwenhoven, D. G. Austing, and S. Tarucha, "Few-electron quantum dots," Rep. Prog. Phys. 64, 701-736 (2001). 
A. Babinski, M. Potemski, S. Raymond, J. Lapointe, and Z. R. Wasilewski, "Fock-Darwin spectrum of a single InAs/GaAs 
quantum dot" Phys. Stat. Sol.(c) 3, 3748-3751 (2006). 

E. Schrodinger, "Der stetige Ubergang von der Mikro- zur Makromechanik" Naturwissenschaften 14, 664-666 (1926). 
E. H. Kennard, "Zur Quantenmechanik einfacher Bewegungstypen" Z. Phys. 44, 326-352 (1927). 

V. V. Dodonov, "Nonclassical states in quantum optics: a squeezed review of the first 75 years" J. Opt.B 4, R1-R33 (2002). 
R. J. Glauber, "Photon correlations" Phys. Rev. Lett. 10, 84-86 (1963). 
U. Leonhardt, Measuring the Quantum State of Light (Cambridge University Press, 1997). 
L. Brown, Quantum Field Theory (Cambridge University Press, 1994). 

D. Leibfried, D. M. Meekhof, B. E. King, C. Monroe, W. M. Itano, and D. J. Wineland, "Experimental determination of 
the motional quantum state of a trapped atom" Phys. Rev. Lett. 77, 281-285 (1996). 
K. Husimi, "Miscellanea in elementary quantum mechanics: I" Prog. Theor. Phys. 9, 238-244 (1953). 
L. Infeld and J. Plebahski, "On a certain class of unitary transformations" Acta Phys. Polon. 14, 41-75 (1955). 

J. Plebahski, "Wave functions of a harmonic oscillator" Phys. Rev. 101, 1825-1826 (1956). 
^® D. Stoler, "Equivalence classes of minimum uncertainty wave-packets" Phys. Rev. D 1, 3217-3219 (1970). 
D. Stoler, "Equivalence classes of minimum uncertainty wave-packets. H" Phys. Rev. D 4, 1925-1926 (1971). 

G. Adesso and G. Chiribella, "Quantum benchmark for teleportation and storage of squeezed states" Phys. Rev. Lett. 100, 
170503 (2008). 

M. Ban, "Quantum dense coding via a two-mode squeezed-vacuum state" J. Opt.B 1, L9-L11 (1999). 
^° S. L. Braunstein and H. J. Kimble, "Dense coding for continuous variables" Phys. Rev. A 61, 042302 (2000). 
S. L. Braunstein, "Squeezing as an irreducible resource" Phys. Rev. A 71, 055801 (2005). 

J. L Cirac, A. S. Parkins, R. Blatt, and P. ZoUer, "'Dark' squeezed states of the motion of a trapped ion" Phys. Rev. Lett. 
70, 556-559 (1993). 

D. J. Heinzen and D. J. Wineland, "Quantum-limited cooling and detection of radio-frequency oscillations by laser-cooled 
ions" Phys. Rev. A 42, 2977-2994 (1990). 

D. M. Meekhof, C. Monroe, B. E. King, W. M. Itano, and D. J. Wineland, "Generation of nonclassical motional states of a 
trapped atom" Phys. Rev. Lett. 76, 1776-1779 (1996). 

E. S. Polzik, "The squeeze goes on" Nature 453, 45-46 (2008). 

^® H. Vahlbruch, S. Chelkowski, K. Danzmann, and R. Schnabel, "Quantum engineering of squeezed states for quantum 
communication and metrology" New Journal of Physics 9, 371 (2007). 

H. Zeng and F. Lin, "Coherent squeezed states of motion in an ion trap generated by Raman-driven sideband transitions" 
Phys. Rev. A 52, 809-814 (1995). 

^® H. J. Kimble, "The quantum internet" Nature 453, 1023-1030 (2008). 

V. V. Dodonov, I. A. Malkin, and V. I. Man'ko, "Coherent states of a charged particle in a time-dependent uniform 
electromagnetic field of a plane current" Physica 59, 241-256 (1972). 

L A. Malkin, V. L Man'ko, and D. A. Trifonov, "Coherent states and transition probabilities in a time-dependent electro- 
magnetic field" Phys. Rev. D 2, 1371-1385 (1970). 

L A. Malkin and V. I. Man'ko, "Coherent states of a charged particle in a magnetic field," Sov. Phys. JETP 28, 527-532 
(1969). 

L A. Malkin and V. I. Man'ko, "Coherent states and excitation of n-dimensional non-stationary forced oscillator" Phys. Lett. 
32A, 243-244 (1970). 

L A. Malkin, V. I. Man'ko, and D. A. Trifonov, "Invariants and the evolution of coherent states for a charged particle in a 
time-dependent magnetic field" Phys. Lett. A 30A, 414 (1969). 

H. R. Lewis and W. B. Riesenfeld, "And exact theory of the time-dependent harmonic oscillator and of a charged particle 
in the time-dependent electromagnetic field" J. Math. Phys. 10, 1458-1473 (1969). 

H. Y. Kim and J. H. Weiner, "Gaussian wave-packet dynamics in uniform magnetic and quadratic potential fields" 
Phys. Rev. B 7, 1353-1362 (1973). 

E. P. Wigner, "On the Quantum Correction For Thermodynamic Equihbrium" Phys. Rev. 40, 749-759 (1930). 
R. P. Feynman and A. R. Hibbs, Quantum Mechanics and Path Integrals (McGraw-Hill, 1965). 



29 



11 



Figure Captions 

Figure 1. x,Px section of Wigner function of squeezed-coherent state, as given by equation (pij) at t = 0, evolving 
under tlie Fock-Darwin Hamiltonian. Tiie lengtli of the x, px axis is, respectively, 10^^ m and 10^^^ Kg m s^^ and 
the Wigner function was multiplied by . 
Figure 2. 

Same as[l]at time t — ■k/{2lob)- 
Figure 3. 

Same as[T]at time t — in / {'iuj b) . 
Figure 4. 

Same as[T]at time t = 57r/(4cj/j;). 
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